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ABSTRACT 

An experimental investigation is reported on slanted base ogive cylin- 
ders at zero incidence. The Mach number range is 0.05 to 0.3. All flow 
disturbances associated with wind tunnel supports are eliminiated in this 
investigation by magnetically suspending the wind tunnel models. The sudden 
and drastic changes in the lift, pitching moment, and drag for a slight 
change in base slant angle are reported. Flow visualization with liquid 
crystals and oil is used to observe base flow patterns, which are respon- 
sible for the sudden changes in aerodynamic characteristics. Hysteretic 
effects in base flow pattern changes are present in this investigation and 
are reported. The effects of a wire support attachment on the 0° slanted 
base model is studied. Computational drag and transition location results 
using VSAERO and SANDRAG are presented and compared with experimental 
results. Base pressure measurements over the slanted bases are made with an 
onboard pressure transducer using remote data telemetry. 
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1. INTRODUCTION 


1.1 Previous Research 

Applications of this study are in the automotive and aircraft industries. The aftbody of 
most transport aircraft fuselages and certain car designs, notably hatchback cars, have slant 
bases. Interest in this research is due to the sudden change in base drag for a small change in 
slant angle. The sudden change in drag is attributed to the sudden change in separated flow 
patterns in the near wake. Specifically, the base flow pattern changes from a quasi-symmetric 
separation pattern to a three-dimensional, longitudinal vortex flow. These vortices produce low 
pressures and therefore, high base drag. 

Jansson and Hucho^ discuss early research on this subject. Their paper presents the effects 
of an automobile rear window on vehicle drag. The study concludes that a small change in rear 
window angle (55" to 65"), creates a large drag change and creates changes in the near wake flow 
pattern. Morel, 2,3 isolated the effects in the near wake by studying slanted base ogive cylinders. 
Fig. 1 illustrates the base flow patterns observed by Morel. This investigation involved the 
study of factors that influence vehicles, such as, rounding the upper slant edge, Reynold’s 
number effects, free stream turbulence, and ground effects. Of these factors, Morel concludes 
that only rounding the upper slant edge had important modifications on the drag overshoot. 
The critical slant angle shifts to a higher angle due to rounding. Drag overshoot existed for all 
factors tested and occurred between 45* and 46.5* for the drag runs and 47.25 and 48 for the 
pressure runs. This variation is attributed to differences in the experimental set-up. Fig. 2 
shows a typical drag coefficient versus slant angle curve. 
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Maull 4 verified Morel’s drag results for slant-based bodies of revolution. Furthermore, drag 
measurements were made for an incidence range of -6 to 6 . Maull shows that the incidence at 
which the drag overshoot occurs is a linear function of the slant angle. Maull also observed a 
hysteresis effect in the formation of the base flow patterns when changing the model incidence. 
Hysteresis effects do not show on drag results since the tunnel was stopped to change incidence 
for drag measurements. 

Xia and Bearman 5 investigated the effects of incidence on axisymmetric bodies with slanted 
bases. Their study shows that bysteretic effects exist for the critical incidence for transition 
between base flow patterns. Furthermore, base pressure measurements with the 45' base at zero 
incidence show that the base flow pattern is vortical for all fineness ratios. However, drag 
measurements for this base at zero incidence indicate that the base flow pattern is a function of 
the fineness ratio. According to Xia and Bearman, these discrepancies may be due to slight 
variations in mounting the models and small variations in the approaching flow properties. A 
further explanation is due to drag measurements being made by increasing the incidence with 
the tunnel running. The authors compare estimates of the incidence at which the base flow 
pattern changes with Maull. Fig. 3 shows the comparison. The lack of agreement is 
attributed to differences in the experimental arrangement (mechanical supports) and 
procedures* 

These earlier investigations of slanted bases used various arrangements for supporting the 
models. Morel obtained measurements by suspending the wind tunnel model on thin wires. 
Maull suspends the wind tunnel model by wires attached to airfoil sections on the model and a 
wire from the straight section of the body. Xia and Bearman used two support arms attached 
to the forebody by two tapered rods, and tail wires. Fig. 4 compares these support systems. 

1.2 Objectives of the Current Study 

Research at the NASA Langley Research Center (LaRC) 13 inch Magnetic Suspension and 
Balance System (MSBS) is intended to validate previous slanted base drag results and to extend 


- 2 - 



lift and pitching moment results. Test results will be free of all flow disturbances due to 
mechanical supports. The slanted base ogive cylinder model in this study is aluminum with 
interchangeable bases. The nose of the model consists of a 7.1 inch ogive radius, while the 
cylindrical midbody is 8.323 inches long with a 1.212 inch outer diameter. The aft body consists 
of interchangeable bases with 0", 30', 40", 45', 50", and 60" slants. A magnetic core of ultra-low 
carbon iron, 7.225 inches long by 1 inch in diameter is inserted in the model for suspension 
purposes. Fig. 5 shows the model. 

1.3 NASA Langley Research Center 13 Inch Magnetic Suspension and Balance System 
1.3.1 Wind Tunnel 

The LaRC 13 inch MSBS wind tunnel is an open circuit design shown schematically in Fig. 
6a. Outside air at ambient conditions is drawn into the wind tunnel through a large bellmouth 
inlet. Air flows through a constant diameter duct from the bellmouth to the first set of 
turning vanes. A quick diffuser and settling chamber are downstream of the first turn. The 
settling chamber has several screens to condition the flow for entry into the contraction and test 
section. A transition section connects the contraction and the test section. The test section 
cross-section is a modified octagon 12.56 inches from top to bottom wall and 10.69 inches 
from side wall to side wall. It is 48 inches long and is made of Lexan (brand name 
Polycarbonate) to permit viewing of the model in flight. Fig. 6b shows the inside dimensions of 
the test section cross-section. Another transition section follows the test section. Both 
transition sections have the same geometry to minimize costs. A 2j inch thick steel honeycomb 
follows the second transition section to protect the fan from loose objects. Another set of 
turning vanes follows the honeycomb to direct the flow into the fan. A 206 hp, 6000 rpm motor 
drives the fan and is installed in a nacelle. 6 
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1.3.2 Electromagnets 

Four copper electromagnets are arranged in a u V-type* configuration to provide the lift 
force, pitching moment, side force, and yawing moment. An axial solenoid upstream of the test 
section provides the drag force. Fig. 6c shows the electromagnet configuration. During the tests 
reported here, the electromagnets were powered by a mix of thyratron, thyristor, and rectified 
motor-driven variac power supplies,^ 


1.3.3 Position Sensors 

The position sensors monitor the model’s position and orientation in five degrees of freedom. 
The system uses five linear photodiode arrays shown schematically in Fig. 6d. Each array 
contains 1024 elements with a 0.001 inch spacing. A laser illuminates a sensor with a vertical 
light sheet. A model in suspension blocks part of the light sheet, casting a shadow on the array. 
The spacing of the elements allows the shadow location to be measured to 0.001 inches 
accuracy. 8 

1.3.4 Control System 

The control system stabilizes and controls the position and attitude of the suspended model. 
The feedback control system is implemented digitally using a PDP 11/23+ minicomputer. The 
controller is essentially a digitial simulation of classical analogue MSBS control systems. 
Programming is almost entirely in assembly language to speed execution. Fig. 6e shows a block 
diagram of the digitally controlled MSBS.^ 


1.4 Data Acquisition 

During testing, the data acquisition system measures the wind tunnel parameters, such as, 
the stagnation and static pressure and the stagnation temperature. Further measurements 
include the coil currents. These currents are converted into aerodynamic forces through a 
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calibration discussed in detail in Chapter 2. 

The Hewlett Packard 3497A Data Acquisition/Control Unit, (DACU) reads thermocouple, 
pressure transducer, and coil voltages. An IBM PC XT runs the software that drives the 
DACU. Calibrations in the software convert voltage readings into temperature, pressure, and 
coil current values. A complete scan through all channels of the DACU measures the control 
and bias currents in the four V-type electromagnets and the control current in the drag solenoid. 
Further measurements within a complete run include the wind tunnel parameters. The data 
acquisition software also computes the average of each measurement using values from the 
present run and from four previous runs. At the completion of each run, the data acquisition 
software calculates the local Mach number using isentropic flow assumptions. 

Modifications to the data acquisition software improve the data retrieval process for future 
research. The DACU scans the coil voltages twenty times before converting these values to coil 
currents and printing them to the screen. The previous software printed these values to the 
screen after every scan which slows the data retrieval process. Pressure transducer voltages are 
read every five runs to calculate the local Mach number. The local Mach number is then 
printed to the screen every fifth scan giving the tunnel operator a faster turnaround in local 
Mach number readings. On the twentieth run, the DACU scans all channels and then prints 
temperature, pressure, current, and Mach number values to the screen. This software is three to 
four times faster than the previous version. 
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2 . ELECTROMAGNETIC FORCE AND TORQUE EQUATIONS 


2.1 Development of Force and Torque Equations 

The magnetic force, F, and torque, T, on a ferromagnetic body can be determined by 
regarding the body as a magnetic dipole placed in a magnetic field. The vector equations can be 
written as follows: 


dF = K t (M-v)HdV (2-1) 

dT=K t (MxH)dV (2-2) 

With reasonably uniform fields, the vector equations can be written as follows:* 


F~K t v(M-V)H (2-3) 

T~K t v(Mxfi) M 

Gilliam*® applied these equations to a MSBS model at a given incidence, yaw, and roll angle. 
Writing these equations in matrix form gives, 


* For example, see An Introduction to Electromagnetic Theory by P. C. Clem mow 
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The components of average magnetization for a soft iron core are, 


Hx 

M y = r+jiu; Ey - 


(2-7a) 


(2-7b) 


(2-7c) 


where fi is the magnetic permeability of the material (relative to air). The relative permeability 
for typical materials that are moderately magnetized is ih the order of 5000 to 20000. The 
demagnetization factors, Dj , are related as follows. 


D x +Dy+D z =l (three dimensional) 


(2-8a) 
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Dy+D z =l (two-dimensional) (2-8b) 

The demagnetizing factor, D x will decrease as the core slenderness ratio increases as Fig. 7 
shows. 11 The MSBS core is cylindrical and has a slenderness ratio of 7.225. For the 
permeability range mentioned above, D x ~ 0.024. Since the relative permeability is quite high, 
the assumption that is allowed. This simplifies the average magnetization equations. 



Expressing the average magnetization equations in the body frame gives: 


M b = cj; Mc= d^ (2 ' 10) 


The average magnetization equations in matrix form gives: 
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Fig. 8 shows the body frame coordinates (a,b,c) with the tunnel reference frame (x,y,z). The 
transformation from the tunnel reference frame to the body reference frame is as follows. 
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( 2 - 12 ) 
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To date research at the LaRC 13 inch MSBS involves zero degree yaw, roll, and angle of attack 
measurements. Therefore, the body frame is aligned with the tunnel reference frame. Under 
these conditions the transformation matrices reduce to the identity matrix as expected. Future 
modifications to the MSBS will allow for high angle of attack tests. Therefore, further analysis 
of the force and torque equations will include the transformation matrix for angles of attack 
with the remaining matrices equal to the identity matrix. 
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(2-13) 


The inverse matrix of a orthogonal matrix is simply the matrix transpose. Therefore, 
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(2-14) 


The transformation matrix allows the average magnetization and the applied magnetic fields to 
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be written in the tunnel reference frame. 



Combining equations 2-9, 2-13, and 2-14 with equations 2-3 and 2-4 gives the following result 
for the force and torque vector. 



Side forces, yawing moments, and rolling moments are not present in this research, therefore, 
simplifying the force and torque matrices give the following equations. 
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(2-19) 
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(2-20) 


,=K t V^^ — ^^|^cos 2 ^-sin 2 ^^H x H z +sintfcos^^H x 2 — H z 2 )| 


Further simplification to the zero incidence case gives: 



( 2 - 21 ) 

( 2 - 22 ) 

(2-23) 

(2-24) 


2.2 Field Assumptions 

Several assumptions are now necessary to fully develop the force and torque equations. 
These assumptions are for the magnetic field and magnetic field gradients which depend on the 
type of magnetic suspension system. The LaRC 13 inch MSBS is a V-type system. The 
assumptions for a V-type system are developed by Basmajian, Copeland, and Stevens 12 and are 
based on magnet symmetry. 
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^-^l( I l l + I l 2 + I l 3 +I l 4 ) 

^=p 2 (I 11 + I 12 + I 13 + I ) 4 )-M 3 I d 

^=A>4( I l 1 - I l 2 + I l3- I l 4 )-^5 I d 

Hx = P6(l| 1 + I l 2 + I l 3 + I l 4 ) + / > 7 I d 

H z =-P 8 (l| 1 -I| 2 +I| 3 -I| 4 ) 


(2-25a) 

(2-25b) 

(2-25c) 

(2-25d) 

(2-25e) 


where through are constants. 

The relative permeability will remain high throughout testing due to the ultra-low carbon iron 
core, therefore, the demagnetizating factors, D a , D^, and D c are constant for a constant volume 
magnetic core. Substituting equations 2-23a through 2-23e into equations 2-20, 2-21, and 2-22, 
realizing that K t is a constant, gives the final results. 


F x = + ^D 2 — 

F 2 = ^sIl 2 + ^iJd — + Vp*D 

Ty = — ^lO^D 

where, 


l L = I l 1 +I l 2 +I| 3 +I| 4 


Ip = L — Ii + Ii — 1| 

P l 2 l 3 i 4 


and S x through <5 10 are constant coefficients. 


(2-26) 

(2-27) 

(2-28) 

(2-29a) 

(2-29b) 


2.3 Determination of Force and Torque Coefficients 

To calculate the coefficients, known forces and torques are applied to the suspended model, 
and the resulting coil currents are recorded. Fig. 9a shows the proper set-up for such a 
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calibration. The pulley arrangement provides the known drag force. The weight pans attached 
to the suspended model inside the test section provide the lift force and pitching moment. This 
calibration procedure could not be undertaken due to the unexpected down-time of the MSBS 
for magnet and power supply modifications. Preliminary calibrations only involved axial loads 
on the 0’ base and the 60* base models. Due to the discovery of significant lift forces and 
pitching moments exerted on the slanted base models, more complete calibration data was 
necessary. Fortunately, the weight and location of the interchangeable bases relative to the 
magnetic core provided built-in lift forces and pitching moments. An attempt to provide further 
lift force and pitching moment data was made by positioning brass rings of a known weight on 
the model. Fig. 9b shows the preliminary calibration set-up. The resulting control and bias 
currents in each coil are recorded for each force and torque data point. 

Following previous analysis, the force and torque coefficients are assumed constant. 
Therefore, the product terms of the force and torque equations can be combined to give the 
following equations. 


— — ^ 4^4 

(2-30) 

Fz=<55X 5 +6 6 X 6 — ^7X7+69^8 

(2-31) 

Ty=£ 9 Xg — ^1()X 10 

(2-32) 


The coefficients can be estimated using multiple linear regression in which the sum of the 
squared residuals are minimized. The range of applied forces and torques cover the range of 
expected aerodynamic forces and torques with limited extrapolation of the data. Table 1. 
contains the final force and torque equations. 
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3 . AERODYNAMIC COEFFICIENTS AND REYNOLDS NUMBER CALCULATIONS 


The aerodynamic coefficients and Reynolds numbers are calculated from thermocouple, 
pressure transducer, and coil current data. Specifically, these measurements are: 

stagnation pressure, p^ 
static pressure, Ps 
stagnation temperature, 
lift current, I| +I| 2 +I| 3 +I| 4 
pitching moment current, 1^ — 
drag current, Iq 

3.1 Aerodynamic Coefficients 

The forces and torques exerted on the model are found from the calibration equations. The 
“wind off” lift, pitching moment, and drag currents supporting the suspended model correspond 
to the no aerodynamic force and torque case. These currents may differ slightly from run to 
run. This is due to small variations in the model position between runs. A simple correction 
easily solves this problem. The change in lift, pitching moment, and drag currents from their 
wind off values are calculated for each data run. These increments are added to the zero force 
and torque calibration currents to give a corrected current value. For each data run, the 
corrected currents are substituted into the calibration equations to give the corresponding lift, 
pitching moment, and drag values. The error in not correcting for the wind off currents is small 
for intermediate force and torque values, however, it becomes a factor for low values. 
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The aerodynamic coefficients are calculated as follows: 



The reference area for the coefficients is the maximum cross-sectional area of the model, while 
the length of the 0’ base is the reference length for the pitching moment coefficient. 


3.2 Reynold’s Number Calculations 

The Reynold’s number based on the model diameter is, 


l 



The absolute viscosity, /iy is found by Sutherland’s Law, 


t*v = 


1.458x 10' 6 T s S 

i + 4?4. 


(3-4) 


(3-5) 


3.3 Buoyancy Corrections 

The thickening of the boundary layer along the test section walls decreases the effective test 
section area. This results in a velocity increase due to mass conservation. Velocity and pressure 


gradients are related by the following equation. 


dPs 

"3x 



(3-6) 
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Fig. 10 illustrates this phenomena. The pressure gradient creates an additional force on a model 
in suspension known as a buoyancy drag. A simplistic equation for calculating buoyancy drag is 
the product of the model volume and the pressure gradient. Glauert 13 modified this approach 
considerably by introducing an effective volume into the correction. Specifically, Glauert shows 
the effective volume of the model can be calculated if the velocity distribution along the model 
is known. 


ad buoy =t 7DT / U Y2dS-A, 7ET (3 ~ 7) 

A panel method computer code called VSAERO determines the velocity distribution, jj, along 
the model length as Fig. 11 shows. Further results from VSAERO are presented in Chapter 5. 
The incremental arc length, dS, equals Rd0 for the ogive nose piece and simply equals dx for 
the cylindrical section. Numerical integration of the above equation results in an effective 
volume of 11.54 in 3 . The actual volume of the ogive cylinder model is 11.39 in 3 , which is 1.3% 
lower than the effective volume. In all cases examined by Glauert, the effective model volume 
was greater than the actual model volume. 

The MSBS test section walls are parallel and therefore do not account for the boundary layer 
growth creating a longitudinal pressure gradient. The procedure for determining the pressure 
gradient in the empty tunnel is as follows. 

Two longitudinal pressure tappings are spaced 12 inches apart, centered around the magnetic 
center. The pressure differential between these tappings are found for various Mach numbers. 
Mach numbers range from 0.002 to 0.35 in increments of approximately 0.05 and then back to 
0.002 with the same increment. A least squares polynomial with the pressure differential as the 
dependent variable and the Mach number as the independent variable fits the data. Fig. 12 
shows the set-up. 
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3.4 Blockage Corrections 

Blockage effects are due to the presence of a body in a closed environment which constricts 
the flow. Conservation of mass requires a blockage to increase the flow velocities around the 
body relative to the unconfmed case. This is evident in Fig. 13. The velocity increment 
between confined and unconfined flows is approximately equal along the model length if the 
body is not “large* compared to the the test section. 14 A correction for the static pressure, 
Mach number, and other flow parameters is necessary to account for the constriction. This 
correction is referred to as “solid body blockage.” 

Due to the expansion of the wake downstream of the model, the velocity outside the wake 
will increase. The test section walls act to increase this velocity more than in the unconfined 
case. Mass conservation requires the velocity around the body to increase to account for this 
effect. Again, a correction is necessary and is referred to as “wake blockage”. Specifically, these 
blockage corrections are the velocity increment due to constricted flow normalized with respect 
to the free stream velocity. They are a function of test section and model geometry. 

The ratio of the ogive cylinder cross-sectional area to the MSBS test section area is 
approximately 1%. This ratio is small and blockage effects are minimal. 15 To achieve as 
accurate results as possible, blockage corrections are made. 

3.4.1 Solid Body Blockage in an Octagonal Tunnel 

Several methods exist for the calculation of solid body blockage for three-dimensional models 
in compressible flow depending on the test section geometry. Batchelor suggests the 
following formula for octagonal test sections. 

■ <m) 

where, 

*.-$/*** . . m 
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Batchelor estimates for moderately small fillets to be sufficiently small in comparison to the 
first term in the squared brackets and is therefore neglected. The MSBS test section fillets are of 
the same order in size and are therefore classified as being sufficiently small. Neglecting this 
term results in a solid blockage equation for a rectangular test section with Tq defined as follows. 


i<m - , * (0 ' 0) (5 - 10) 


Carrying the summation from 500 to -500 results in a converged solution of r 0 — 0.82 and a solid 
body blockage factor of e*=7.576x!0 3 


3.4.2 Wake Blockage in an Octagonal Tunnel 

For bodies of revolution in a rectangular tunnel, Allen and Vincenti and Herriot 18 suggest 
the following formula for wake blockage. 



It is further suggested to replace the rectangular area, BH by the octagonal test section area, C. 
3.5 Corrections to Flow Properties 

Corrections to flow properties due to blockage can be made now that the blockage factors, 
and e s exist. As stated, these factors combine to give the total non-dimensional velocity 
increment due to constrictions in a closed-throat tunnel. Mathematically, this gives the velocity 
correction as follows. 

U c = U+AU b (3- 12 ) 

or, ... 
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U c =(l+€ t )U where, e t =e s + f w 


(3-13) 


For compressible flow, the Mach number increment is of importance. Differentiation of the 
isen tropic flow equation, 





(3-14) 


gives the following increment. 


2 dU = dM (3-15) 

^ U 3 M 3 


But, C 


_dU 


t“TT 


by definition. 


Therefore, 


ao 2 « t M 


( 2JlV .M= dM 

\V 2 J \* 2 J 


(3-16) 


Substitution of Equation 3-14 into Equation 3-16 gives the final result. 




(3-17) 


The corrected Mach number becomes, 

M c =M+(l+^pM 2 )Me t (3-18) 

A similiar result exists for the static pressure increment upon differentiation of the isentropic 
flow equation. 
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(3-19) 


Pt 

Pi 




The static pressure increment is as follows. 


dp8 = -7M 2 Pse t 


giving, 


Ps c =P8-7M 2 p 8 c t 


Similarly, the corrected static temperature is as follows. 




(3-20) 


( 3 - 21 ) 


(3-22) 


Fig. 14 compares buoyancy and blockage drag corrections with uncorrected values. Buoyancy 
corrections are between 1.5% and 6% of the uncorrected drag and are the largest corrections to 
the drag data. Blockage corrections are between 1% and 2% of the uncorrected drag. 
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4 . EXPERIMENTAL RESULTS 


4.1 Aerodynamic Characteristics 

Experimental results from the LaRC MSBS are presented in Figs. 15a through 15r. MSBS 
results are corrected for buoyancy and blockage by the methods presented in the previous 
chapter. Scatter exists at the lower Reynolds numbers due to the inability to accurately resolve 
current and pressure readings at very low speeds. This is evident in most of the results 
presented. The reference area for the drag, lift, and pitching moment coefficients is the 
maximum cross-sectional area of the model. The reference length of the pitching moment 
coefficient is the length of the 0* base model. 

Several runs with the 0* base exhibit fair repeatability of drag results. Lift and pitching 
moment data are not well resolved due in part to the imperfect calibration data available. 
Available computational results presented in Chapter 6. and flow visualization results discussed 
at the end of this chapter indicate the boundary layer is laminar over the entire body for these 
lower Reynolds numbers. It is therefore reasonable that runs with transition fixed just 
downstream of the nose piece indicate a substantial rise in drag due to the higher stresses 
associated with a turbulent boundary layer. 

Free transition on the 30* and 40* bases show that lift and pitching moment results increase 
proportionally with drag. This is thought to be evidence of low pressures on the base as 
illustrated in the sketch below. 


Lift force 


Resultant force 


Drag force 



pitching moment 


- 22 - 



An unexpected drag increase occurred with the 45 base with increasing Reynolds number. 
This increase still exists upon lowering the Reynolds number past its critical value implying a 
hysteretic effect. A similiar lift and pitching moment increase is coincident with drag results. 
The sudden change in aerodynamic forces and torques appear to be due to sudden changes in 
base flow patterns. A repeat run was made to more closely approximate the critical 
Reynolds number where the increase occurred. The critical Reynolds number (Rcq) is 
approximately 60,000. An attempt to explore the effects of fixing transition on the 45* base 
was then made with transition fixed just downstream of the nose piece. It was found that the 
sudden drag increase did not occur for Reynolds numbers greater than the critical value with 
fixed transition. Similiar results apply to the lift and pitching moment. Previous research by 
Page 18 on axisymmetric body base flow patterns concludes that base flow patterns are a strong 
function of the incoming boundary layer. Differences in the base flow pattern exist on the 45 
base due to fixed and free transition. This implies that this conclusion may be extended to non- 
symmetric bases. 

Data for the 50* and 60* bases could not be taken at higher speeds due to some stability 
problems (uncontrolled rolling of the 60* model) and large vertical loads on the rear of the 
model (reducing the rear electromagnet currents). 

Fig. 16 illustrates a comparison of drag results between Morel and MSBS results interpolated 
from experimental data for Re D =94000. Morel’s results for the lower base slant models lie 
between fixed and free transition MSBS results. The wire supports used by Morel are likely to 
have caused partial turbulent boundary layers responsible for these drag discrepancies. Existing 
measurements and computational results indicate the boundary layer is laminar over the entire 
body for Re o =94000. Morel measured the drag overshoot at a slightly higher slant angle of 
46.5* (for the drag runs), while MSBS data indicate that the drag overshoot is Reynolds number 
dependent and exists for the 45* base. Time limitations prevented the testing of a cluster of 
bases around 45* to further resolve the overshoots. There is good agreement between MSBS 
data and results by Morel for the 50* and the 60* bases. Fig. 17 and Fig. 18 show the 
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corresponding lift and pitching moment results versus base slant angle for Re D =94000 
respectively. 

4.2 Flow Visualization 

The effects of a wire attachment on an ogive cylinder model was made. A piece of grit 
(approximately 0.5 mm in diameter) was placed just downstream of the nose piece to represent 
the wire attachment. Flow visualization liquid crystals were applied to the midbody. Liquid 
crystals selectively reflect discrete wavelengths (color) of light in response to a shear stress.^® 
Liquid crystals can therefore distinguish between laminar and turbulent boundary layers. Fig. 
19 shows flow visualization results for the modeled wire attachment. The spotted areas on the 
nose piece are color changes representing areas of high stresses. A region of high stresses also 
develops downstream of the grit representing a turbulent wedge. This result further explains 
why Morel’s drag results lie between fixed and free transition MSBS results. 

An attempt at flow visualization using liquid crystals and oil flow was made over the slanted 
bases. Fig. 20a shows the quasi-symmetric flow pattern over the 30' base and shows pressure 
distributions by Xia and Bearman. 5 A very distinct horseshoe pattern exists over the 50’ base 
representing the vortical flow pattern. Fig. 20b illustrates this pattern along with pressure 
distributions. The line surrounding the horseshoe corresponds to the region of minimum 
pressure, or the vortex center. The darker colored region inside this line appears to be a region 
of reattachment. Oil flow over the 45' base revealed the change in base flow patterns at the 
critical Reynolds number. The flow pattern remains vortical for decreasing Reynolds numbers 
past its critical value confirming the hysteretic effect. 

Further attempts at flow visualization using liquid crystals over the cylindrical midbody was 
made for fixed transition just downstream of the nose piece. The characteristic liquid crystal 
color change occurred, representing a high stress region (turbulent boundary layer) not present 
on the clean midbody. 
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5 . COMPUTATIONAL RESULTS 


Two codes are used for predicting transition locations and drag for the (T slanted base 
model, SANDRAG 20 , predicts drag of bodies of revolution at zero angle of attack in 
incompressible flow. VSAERO 21 , a surface panel method for predicting subsonic aerodynamic 
characteristics of arbitrary configurations is also used. 

5.1 Transition 

Fig. 21 shows predicted boundary layer transition and surface pressure distributions for the 
0* slanted base. The transition locations predicted by VSAERO and SANDRAG are in poor 
agreement. VSAERO predicts transition far forward of SANDRAG predictions. Present 
experimental results using flow visualization agree with SANDRAG results in that transition 
starts at the back end and moves forward with increasing Reynolds number, 

To understand the reasons for these results, the factors that influence transition need to be 
analyzed. VSAERO predicts transition in the region of adverse pressure gradient for the 
Reynolds numbers tested. In this region Tollmien-Schlichting instability grows rapidly. In this 
case VSAERO predicts the instability is sufficiently large to cause transition. It is of interest 
that a large portion of the body has a flat pressure gradient. This may cause a rapid growth in 
the Tollmien-Schlichting instability creating a considerable amount of uncertainity in the 
prediction of transition location in this region. VSAERO can predict accurately the transition 
location of bodies with strong favorable pressure gradients in which Tollmien-Schlichting waves 
are dampened. 

SANDRAG relies on built-in empirical correlations to predict transition on axisymmetric 
bodies. This code is able to consider surface roughness in its calculations. Computed transition 
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locations on the configuration vary considerably with the surface roughness. 

5.2 Drag Results 

Most of the drag on this configuration is due to base pressure drag. SANDRAG predictions 
use empirical results to determine the base pressure drag component. VSAERO is able to 
determine base pressure drag by attaching a wake to the base. Base pressures computed by 
VSAERO are a strong function of the tangential velocities inside and outside the vortex sheet 
wake. There is considerable difficulty in choosing a good approximation of these parameters and 
for this reason, a drag comparison between VSAERO and SANDRAG is excluded. 

Fig. 22 shows drag variations with Reynolds numbers using SANDRAG. Predictions 
indicate that a majority of the drag on this configuration is due to base pressure drag which is 
presumed constant (C D BASE=0.13) for all Reynolds numbers tested. The drag increases for 
Re D greater than 2.6xl0 5 due to transition moving forward from the model base. Flow 
visualization results show that transition does occur at the model base and moves forward with 
increasing Reynolds number. Table 2 compares experimental and computational results for 
Re D =94000. 

Fig. 23 compares SANDRAG predictions with MSBS results for the 0* base. There is good 
agreement between MSBS wind tunnel results and SANDRAG predictions for the lower 
Reynolds numbers, however, MSBS results indicate that transition occurs sooner than 
SANDRAG predictions. Differences in transition locations can be attributed to variations in 
surface roughness and differences in turbulence intensities between MSBS wind tunnel values and 
the values of the empirical correlations used in SANDRAG. The slanted base models have a 
standard machine surface roughness of 10 /iin which lies between the two surface roughness cases 
of Fig. 22. The 13* MSBS wind tunnel has a turbulence intensity < 0.1%. The Appendix 
compares calculation methods between SANDRAG and VSAERO. 
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Table 2 0' base drag comparison, Re D =94000 


Morel 

MSBS 

SANDRAG 

■ 

clean 

fixed 

0.186 

m 

n 
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6 . BASE PRESSURE MEASUREMENTS 


MSBS’s eliminate the need for mechanical supports, therefore, another method of retrieving 
base pressure measurements is necessary. Wireless data acquisition using an infrared pressure 
telemetry system solves this problem. The telemetry system contains a pressure transducer, a 
signal conditioning circut, an infrared LED (Light Emitting Diode), and a power source. This 
system was developed for base pressure measurements by Tcheng, Schott, and Bryant 23 , and is 
installed in the wind tunnel model. 

6.1 Pressure Telemetry System 

Fig. 24 shows a schematic of the pressure telemetry system. The pressure transducer is an 
Endevco 8510B 2 psid piezoresistive differential pressure transducer. The vent tube of the 
transducer is filled with argon gas at atmospheric pressure and is sealed with epoxy. This is a 
reference pressure. The transducer contains a four-arm strain gage bridge diffused in a silicon 
diaphram for maximum sensitivity. Fig. 25 shows the dimensions of the pressure transducer. 

The signal conditioner converts the analog output of the pressure transducer into a pulsed 
optical signal. It contains a LM124 quad operational amplifier for signal conditioning. Three of 
the operational amplifiers differentially amplify the transducer signal which operates on very low 
power. An AD537 voltage to frequency converter is used to produce a square wave output from 
its DC voltage input. The input stage of the AD537 converts the input voltage into a drive 
current. The drive current charges a timing capacitor which controls the astable multivibrator 
frequency. The square wave output of the AD537 drives a NPN transistor (2N2222) which 
charges another capacitor. This capacitor discharges through the LED (OP160W) when the 
square wave is high. The LED has a flat window providing a wide radiation angle. This allows 


- 28 - 


for ease in data transmission. The pulsed optical signal from the LED is received by a silicon 
photodiode located outside the test section and the frequency is read by a digital frequency 
counter. 


The power supplies for the pressure telemetry system are hearing-aid batteries (E41E) which 
are rated at 1.4V each. Five batteries in series provide the power requirements. 

6.2 Wind Tunnel Pressure Model 

A slanted base pressure model was developed to house the telemetry package. The major 
problem in installing the pressure telemetry system in the model is space restriction. 


6.2.1 Installation of the Pressure Telemetry System 

The signal conditioner flies in the nose piece of the model due to size limitations of its 
components. It is housed in a 0.75 inch diameter hole drilled to a depth of 1.6 inches. The 
signal conditioner is built on both sides of a circut board in order to minimize its length. The 
battery pack flies forward of the signal conditioner in a 0.5 inch diameter hole drilled to a depth 
of 0.96 inches. Fig. 26a shows the model nose piece for pressure measurements. 

The pressure transducer and the LED are housed in the model aftbody. A 0.125 inch 
diameter hole is drilled through the magnetic core to run connections between the signal 
conditioner, the pressure transducer, and the LED. Fig. 26b shows the magnetic core used for 
pressure measurements. A 0.17 inch diameter hole through the side of the aftbody houses the 
LED. The LED is epoxied in place with its window flush against the outside of the aftbody. 
The pressure transducer has a 10-32 mounting thread which screws into the rear of the aftbody. 
The face of the transducer is flush with the rear of the aftbody after installation. Fig. 26c shows 
the aftbody which houses the pressure transducer and the LED. The slanted bases fit onto a 
shaft located on the rear of the aftbody. A RC1 fit is specified to join the aftbody and the 
slanted bases. Rubber cement, is applied at the joint to ensure an air-tight fit. Fig. 26d shows 
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the slanted base design for pressure measurements. Notice that a cylindrical pocket 1.131 inches 
in diameter and 0.025 inches long exists between the aftbody and the slanted bases. The need 
for this pocket will be discussed later. 

Pressure taps are drilled on the 0*, 40*, 45', and 50* bases. Flow visualization results from 
Chapter 4. show the longitudinal vortex is centered near a distance of | the radius (r) from the 
base center. A horizontal pressure tap pattern on the slanted bases is chosen with a geometric 
spacing of Or, |r, |r, and |r to show the horizontal pressure gradient. A similiar pressure tap 
pattern is chosen to measure the vertical pressure gradient. Fig. 26d also shows the pressure tap 
pattern. 

6.2.2 How the Pressure Telemetry System Works 

To measure the pressure at a specific station, all pressure stations not in use are sealed. The 
pocket between the aftbody and the slanted base is open to the outside by the unsealed pressure 
tap. At steady state, the pocket will reach the pressure at the unsealed station. The pressure 
transducer actually measures the pressure in this pocket that corresponds to this station. 

6.3 Calibration of the Pressure Transducer 

Fig. 27 shows a schematic of the pressure calibration set-up. The calibration procedure 
involves applying known pressures to the pressure transducer at constant temperature and 
measuring the resulting frequency output across the LED. The pressure is varied by a pressure 
regulator attached to the pressure transducer. A manometer is in line with the pressure 
regulator to measure the resulting pressures. The LED signal is measured by a photodiode in 
line with a frequency counter. Gage pressures are varied between 0.0 and 20.0 in. I^O and the 
resulting LED signal is recorded. Fig. 28 shows that pressure is a linear function of the square 
wave frequency at constant temperature. Fig. 28 also shows the pressure transducer is extremely 
sensitive to temperature changes. The pressure transducer is in good thermal contact with the 
aluminum model shell, and the transducer face is exposed to tunnel conditions through the 
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pressure taps. The approximation is therefore made that the pressure transducer temperature is 
equal to the tunnel stagnation temperature. 

Fig. 29 shows the effects of varying the input voltage on the frequency output for two 
limiting temperatures. The telemetry system does not have a voltage regulator, and therefore, 
the battery voltage will decrease with operation. Battery voltages were measured before and 
after operation to determine the operational range. The output frequency is shown to be 
independent of the input voltage over the entire operational range of the batteries. 


6.4 Experimental Results 

Base pressure coefficients are determined by the following equations: 

^{P(base)“P 8 } ft i\ 

Cp(base) = 7y v 6 * 1 ) 

P TPsM 2 

where, P(base) = p t( M=0 ) Ptransdu cer 

Ptransducer ™ f° und froin the calibration curves for a given output frequency and transducer 
temperature. Pressure measurements were made on the 40* and 50* base along a horizontal line 
passing through the centerline. Results for the 40* and 50 base are shown in Fig. 30. Scatter 
exists in MSBS measurements due to the sensitivity of the pressure transducer to temperature 
changes. 

There is good agreement between MSBS (ReQ= 100,000) and Morel’s 24 results 
(Re D =94,000) on the 40* base for the centerline pressure station. Xia and Bearman’s 5 results 
(Re o =180, 000-290, 000) consistently lie above MSBS results implying there may be a Reynolds 
number effect on base pressure. The pressure gradient trends between MSBS and Xia and 
Bearman’s results are similiar, showing the flow pattern for the 40 base is quasi-symmetric. 
Morel only documents the centerline pressure value for the 40 base, therefore, a pressure 
gradient cannot be deduced. 
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There is good agreement between MSBS and Morel’s results on the 50 base. As with the 
40* base, Xia and Bearman’s base pressure measurements consistently lie above MSBS results 
implying a Reynolds number effect. Large pressure gradients exist for all results, showing the 
vortical flow pattern is present on the 50* base. MSBS results show the pressure maxima occurs 
near 2^ = |. This agrees with MSBS flow visualization results and Xia and Bearman’s results. 
Morel only documents three pressure values, therefore, a horizontal pressure maxima cannot be 
accurately deduced. 

A source of error in MSBS base pressure measurements is due to the inability to accurately 
resolve the pressure transducer temperature during operation. The transducer temperature is 
assumed to be equal to the wind tunnel stagnation temperature during operation. When not in 
use, the base pressure model was stored at room temperature, which was 5 to 10 F below the 
wind tunnel stagnation temperature. The pressure model was suspended in the wind tunnel at 
low Reynolds numbers for 10-15 minutes before measurements were taken in order to stabilize 
the pressure transducer temperature. Thermal effects of the iron core and the aluminum model 
may still be present after this time, causing the pressure transducer temperature to be less than 
the stagnation temperature. Electrical heating may also be present during operation which will 
tend to increase the transducer temperature. The magnitudes of these thermal effects are not 
known. Inaccuracies in the transducer temperature measurement can be minimized by direct 
measurement with a thermocouple, however, this creates the need for temperature telemetry. 

Further problems develop in determining the pressure differential, P(base) — Ps ‘ This 
quantity is not measured directly, therefore, the pressure differential is found by subtracting two 
large numbers. The pressure differential is between .1% and 2% of these pressure values. 
Solving this problem involves referencing the wind tunnel stagnation pressure instead of the 
pressure of the sealed argon gas. This will also eliminate the temperature assumption necessary 
in this study. 


- 32 - 



7 . CONCLUSIONS 


The sudden change from a quasi-symmetric base flow pattern to a vortical base flow pattern 
has been shown to be responsible for the large drag increase on slanted base ogive cylinders at 
zero incidence. Flow visualization using liquid crystals and oil flow shows the base flow pattern 
change which further validates MSBS drag results. Lift forces and pitching moments acting on 
the slanted base configurations were also measured revealing large increases similiar to the drag 
results. 

Previous research^ shows that the change in base flow patterns is not “too dependent on 
Reynolds number, however, this research shows the base flow pattern change occurs suddenly 
with increasing Reynolds number (Rep ~ 60, 000) on the 45 base. Hysteresis is present in base 
flow pattern changes on the 45’ base with decreasing Reynolds number. Tripping the boundary 
layer far forward on the 45* base prevents the formation of the vortical flow pattern and 
therefore prevents any hysteretic effects. A summary of MSBS drag results with Morel’s drag 
results shows that Morel’s results lie between MSBS fixed and free transition results. Flow 
visualization with liquid crystals reveals that partial turbulent boundary layers form due to 
wire support attachments on the model. 

Computational results are difficult to achieve for the 0* slanted base configuration unless 
empirical results are used. This is primarily due to the large contribution of base drag, 
responsible for more than 50% of the total drag, which is not easily modeled computationally. 
Transition locations are not easily modeled computationally due to the flat pressure gradient 
over most of this configuration. In order to achieve computational results for the slanted base 
configurations a three-dimensional Navier Stokes solution is needed. 

Base pressure measurements on the 40* and 50* bases using remote data telemetry show the 
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two flow patterns present in this research. In conclusion, Magnetic Suspension and Balance 
Systems can retrieve interference-free data on aerodynamically complex configurations not 
available in conventional wind tunnels. 
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APPENDIX A: FUTURE RESEARCH 


Aerodynamic Testing 

Previous research at the LaRC 13 inch MSBS by Dress 25 showed the ability of the system to 
take force measurements. This research involved drag measurements on bodies of revolution at 
zero incidence. No lift forces or pitching moments occurred for this test case, therefore, no 
interactions between coil currents developed. The drag force is a function of the current in the 
drag solenoid, I D . 

Large lift forces and pitching moments were shown to act on slanted base ogive cylinders at 
zero incidence creating significant coil current interactions. The development of the zero 
incidence equations for drag, lift, and pitching moment in Chapter 2. account for these 
interactions. These equations allow for complex configurations to be tested at zero incidence. 

The next step in MSBS research will be the aerodynamic testing of complex configurations 
at arbitrary angles of attack. Lift, pitching moment, and drag measurements for a fighter 
configuration at a given angle of attack is an example of such research. The governing 
equations for this research are determined by substituting the field assumptions (2-23a through 
2-23e) into the angle of attack equations (2-17 through 2-19). The coefficients are constant for a 
constant angle of attack and can be determined by the calibration precedure outlined in Chapter 
2 . 


Pressure Measurements 

Chapter 6. shows how pressure measurements are possible for magnetically suspended models 
using remote data telemetry. The calibration shows that the pressure transducer is extremely 
sensitive to temperature. An onboard thermocouple with temperature telemetry is necessary to 
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accurately determine the pressure transducer temperature and will eliminate the temperature 
approximation necessary in this research. Furthermore, the pressure differential, P(base)'P 8 ’ 
needs to be measured directly for accurate results. Some barriers still exist preventing multiple 
pressure measurements. This research has shown how multiple pressure measurements can be 
made over a small surface area (in this case the slanted bases) with a single channel pressure 
transducer. However, before multiple pressure measurements over large surface areas can be 
made, a miniature multi-channel pressure transducer needs to be developed. 
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APPENDIX B: CALCULATION METHODS IN SANDRAG AND VSAERO 


The two codes used in this research, VSAERO and SANDRAG, implement various methods 
to solve for the flow fields, SANDRAG results are restricted to axisymmetric bodies at zero 
incidence in incompressible flow, however, this code is “user-friendly". VSAERO can predict 
subsonic flow fields for general configurations at any orientation, however, this code is not user- 
friendly. 

The effects of viscosity in VSAERO are determined in an iterative loop which couples the 
potential flow and boundary layer solutions. VSAERO is a panel method code in which doublet 
and source singularities are piecewise constant over the panels,^* Thwaites method with 
Curie’s modifications is used in VSAERO for laminar boundary layer calculations, while 
Granville’s criterion is used for transition predictions. Nash and Hicks method is used for 
turbulent boundary layer calculations. Wake calculations in VSAERO determine base pressure 

• 22 
estimates. 

As with VSAERO, viscous effects are determined by coupling the potential flow and 
boundary layer solutions. SANDRAG uses an axial distribution of source and sink elements to 
determine the potential flow solution. These strengths vary linearly with length. Laminar 
boundary layer solutions in SANDRAG are determined by Thwaites’ one parameter solution of 
the momentum integral equation for steady flow over axisymmetric bodies. Transition 
calculations in SANDRAG are determined by an approximate method by Schlichting and 
Ulrich, and by combining surface roughness effects using Kluck’s empirical results. Turbulent 
boundary layer calculations are made by a modification of White’s Karman-type method which 
includes surface roughness effects. Base pressure estimates in SANDRAG are found by empirical 
results by Payne, Hartley, and Taylor for various axisymmetric afterbodies . 20 
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(a) Qua si-axisym metric Separation Pattern 
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Fig. 3 Critical angle of attack for drag overshoot. 




(b) Maull 



Fig. 4 Support systems (not to scale). 
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Fig. 5 Slanted base ogive cylinder model components. 







(c) Electromagnet configuration 



(d) Position sensors 


Fig. 6 cont’d 




(e) Block diagram of 13" MSBS hardware 


Fig, 6 cont'd 
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Fig. 7 Demagnetizing factor for a cylinder, (Dx) . 





(a) Proper calibration set-up 
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(b) Flow with boundaries. 

Fig* 10 Buoyancy effects in a closed throat wind tunnel* 
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Magnetic 
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Fig. 12 Pressure tap locations to determine pressure differential. 




(a) Unconfined flow 



Fig. 13 Blockage effects in a closed throat wind tunnel. 




Uncorrected 
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Fig. 14 Buoyancy and blockage corrections on drag results. 

(0° base) 
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(b) Lift and pitching moment, 0° base, 

Run 001/003, free transition 

Fig. 15 MSBS drag, lift, and pitching moment results. 
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Fig. 15 cont'd 
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Fig. 15 cont'd 
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Fig. 16 MSBS/Morel drag results. 
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Fig. 18 MSBS pitching moment results. 
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Fig. 19 Effects of a wire attachment. 








2. Pressure distribution (50° base). 

Fig. 20b Flow visualization on 50° base (Liquid crystals). 
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Transition (VSAERO) Transition (SANDRAG) 
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Fig. 21 Pressure distribution (VSAERO) and predicted 
transition locations (0° base) . 



Fig. 22 Details of SANDRAG computations. 




Fig. 23 Comparison of SANDRAG and MSBS drag results. 
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26 Details of pressure telemetry mode 1 . 
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(c) Aftbody 
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(d) Slanted bases 
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Fig. 30 MSBS slanted base 
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